Casimir scaling in gauge theories with a gap. Deformed QCD as a toy model. 
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We study a Casimir- like behaviour in a "deformed QCD". We demonstrate that for the system 
defined on a manifold of size L the difference A_B = E — -Biviink between the energies of a system 
in a non-trivial background and Minkowski space-time geometry exhibits tire Casimir-like scaling 
AE ~ L~^, despite the presence of a mass gap in the system, in contrast with naive expectation 
AE ~ exp(— mL), which would normally originate from any physical massive propagating degrees 
of freedom consequent to conventional dispersion relations. The Casimir-like behaviour in our 
system comes instead from a non-dispersive ("contact") term which is not related to any physical 
propagating degrees of freedom, such that the naive argument is simply not applicable. These ideas 
can be explicitly tested in weakly coupled deformed QCD. We comment on profound consequences 
for cosmology of this effect if it persists in real strongly coupled QCD. 



I. INTRODUCTION AND MOTIVATION 

The main motivation for the present studies is a recent 
suggestion on the dynamical Dark Energy (DE) model 
which is entirely rooted in the strongly coupled QCD, 
without any new fields and/or coupling constants [1- 
4]. The key element of the proposal [1-4] is based on 
paradigm that the relevant energy which enters the Ein- 
stein equations is in fact the difference AE = E — E'lviink 
between the energies of a system in a non-trivial back- 
ground and Minkowski space-time geometry, similar to 
the well known Casimir effect when the observed energy 
is a difference between the energy computed for a system 
with conducting boundaries (positioned at finite distance 
L) and infinite Minkowski space ^. This paradigm is 
based on the conjecture that gravity, as described by the 
Einstein equations, is a low-energy effective interaction 
which, as such, should not be sensitive to the microscopic 
degrees of freedom in the system but to some effective 
scale. Thus, the energy density that enters the semiclas- 
sical Einstein equations should not be the "bare" energy 
as computed in QFT, but rather a "renormalized" energy 
density. We propose the renormalization scheme given 
above which sets the vacuum energy to zero in Minkowski 
space wherein the Einstein equations are automatically 
satisfied as the Ricci tensor identically vanishes. 

The above prescription is in fact the standard sub- 
traction procedure that is normally used for the descrip- 
tion of horizon thermodynamics [5, 6] as well as in a 
course of computations of different Green's function in 
a curved background by subtracting infinities originated 
from the flat space [7]. In the present context such 
a definition AE = {E — i^Mink) for the vacuum en- 
ergy was first advocated in 1967 by Zeldovich [8] who 



^ Here and in what follows wc use term "Casimir effect" to em- 
phasize the power like sensitivity to large distances irrespectively 
to their nature. A crucial distinct feature which characterizes 
the system we are interested in is the presence of dimensional 
parameter L ~ H~^ (where H is a Hubble constant) in the 
system which discriminates it from infinitely large Minkowski 
space-time. 



argued that pvac ~ Gm2, with rUp being the proton's 
mass. Later on such a definition for the relevant en- 
ergy AE = {E — Euink) which enters the Einstein equa- 
tions has been advocated from different perspectives in a 
number of papers, see, for example, the relatively recent 
works [9-16] and references therein. 

We study the scaling behavior of AE when the back- 
ground deviates slightly from minkowski space. The 
difference AE must obviously vanish when any devia- 
tions (parametrized by Hubble constant or inverse size 
of the visible universe, iJ ~ L^^) go to zero as this cor- 
responds to the transition to flat Minkowski space. A 
naive expectation based on common sense suggests that 
AE ~ exp(— Aqcd/-H^) ^ exp(— 10'*^) as QCD has a mass 
gap ^ Aqcd ~ 100 MeV, and therefore, AE must not be 
sensitive to size of our universe L ^ H~^. Such a naive 
expectation formally follows from the dispersion relations 
which dictate that a sensitivity to very large distances 
must be exponentially suppressed when the mass gap is 
present in the system^. 

However, as emphasized in [3, 4] in strongly coupled 
gauge theories along with conventional dispersive con- 
tribution there exists a non-dispersive contribution, not 
related to any physical propagating degrees of freedom. 
This non-dispersive (contact) term generally emerges as 
a result of topologically nontrivial sectors in four dimen- 
sional QCD. The variation of this contact term with vari- 
ation of the background may lead to a power like scaling 
AE ^ H + 0{H)^ rather than to an exponential like 
AE ~ exTp{— Aqcu / H) . If true, the difference between 
two metrics (FLRW and Minkowski) would lead to an 
estimate 



AE 



A3 



(lO^-'cl/)'*, 1/L - i? - IQ-^'^eF (1) 



which is amazingly close to the observed DE value today. 
It is interesting to note that expression (1) reduces to 



^ The Casimir effect due to the massless EkM field obvious shows 
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such power dependence A_E = — JI 4 . Similar computations for 
a massive scalar particle with mass m leads to an exponentially 
suppressed result AE ~ exp(— mL) as expected, see e.g. [17]. 



Zeldovich's formula pvac ^ Grrip if one replaces Aqcd ^■ 
nip and H — >■ GAq^^. The last step follows from the 
solution of the Friedman equation 



2 SttG , , 

J^ = -^r- [PDE + Pm) , Pde 



^Aqcd 



(2) 



when the DE component dominates the matter compo- 
nent, PDE ^ PM- In this case the evolution of the uni- 
verse approaches a de-Sitter state with constant expan- 
sion rate H ~ GAq^jj^ as follows from (2). 

Another motivation to study the Casimir like be- 
haviour in QCD is a proposal [18, 19] that the V odd 
correlations observed at RHIC and LHC is in fact another 
manifestation of long range order advocated in this work. 
Furthermore, an apparently universal thermal spectrum 
observed in all high energy collisions when the statisti- 
cal thcrmalization could never be reached in the systems, 
might be also related to the same contact term, not re- 
lated to any physical propagating degrees of freedom, see 
[18, 19] and references therein for the details. 

There are a number of arguments supporting the power 
like behaviour AE ^ H + 0{H)^ in gauge theories, see 
section IV where we present some general arguments sug- 
gesting the Casimir like corrections in gauge theories with 
nontrivial topological structure. However, it is always de- 
sirable and very instructive to sec how the general argu- 
ments work in some simplified settings. This is precisely 
the goal for the present study: we want to consider a sim- 
plified ("deformed") version of QCD which, on one hand, 
would be a weakly coupled gauge theory wherein com- 
putations can be performed in theoretically controllable 
manner. On other hand, this deformation would preserve 
all the relevant elements of strongly coupled QCD such as 
confinement, degeneracy of topological sectors, nontrivial 
6 dependence, presence of non-dispersive contribution to 
topological susceptibility, and other crucial aspects, for 
this phenomenon to emerge. Such a "deformed" theory 
has recently been developed [20]. All computations in 
this work (excluding those in sec. Ill A) are performed 
within this framework. 



II. TOPOLOGICAL SUSCEPTIBILITY IN THE 
DEFORMED QCD 

In the deformed theory an extra term is put into the 
Lagrangian in order to prevent the center symmetry 
breaking that characterizes the QCD phase transition 
between "confined" hadronic matter and "deconfined" 
quark-gluon plasma. Thus we have a theory which re- 
mains confined at high temperature in a weak coupling 
regime, and for which it is claimed [20] that there does 
not exist an order parameter to differentiate the low tem- 
perature confined regime from the high temperature con- 
fined regime. The non-trivial topological sectors of the 
theory are described in this model in terms of the weakly 
coupled monopolcs. The monopoles in this framework 
arc not real particles, they are pseudo-particles which live 



in Euclidean space and describe the physical tunnelling 
processes between different topological sectors \n) and 
|n-|- 1). In particular, the monopole fugacity C. should be 
understood as number of tunnelling events per unit time 
per unit volume 



/number of tunnelling events \ NcC 
[ VL J " "T' 



(3) 



where extra factor Nc in (3) accounts for Nc different 
types of monopoles present in the system and L is the 
size of the circle along t = it and plays the role of the in- 
verse temperature. The monopole gas experiences Debye 
screening so that the field due to any static charge falls off 
exponentially with characteristic length m~^. The num- 
ber density Af of monopoles is given by the monopole 
fugacity, ^ C, so that the average number of monopoles 
in a "Debye volume" is given by 



M = 



x-^ 



1 



»1, 



(4) 



The last inequality holds since the monopole fugacity is 
exponentially suppressed, ( ^ e^^'^ , and in fact we can 
view (4) as a constraint on the validity of the approxi- 
mation. 

The topological susceptibility in this model can be ex- 
plicitly computed and is given by [21] 



XYM = / d^x{q{x)q{0)) = 



C 



NrL 



d^x [6{k)] . (5) 



The light quarks can be easily inserted into the system. 
The corresponding generalization of eq. (5) reads [21] 



XQCD 



d^x{q{^)q{0)) 



(6) 



C 



NrL 



d^s 



5{X.) - 7X1%, 



Anr 



The first term in this equation has non-dispersive nature, 
similar to eq. (5) and has the positive sign. This contact 
term (which is not related to any physical propagating 
degrees of freedom) has been computed in this model 
using monopoles describing the transitions between the 
degenerate topological sectors '^. The positive sign of 
this term is the crucial element of the resolution of the 
U{1)a problem. The second term in eq. (6) is a standard 



^ In the context of this paper the "degeneracy" inipUes there exis- 
tence of winding states \n) constructed as follows: T\n) = |n-|-l). 
In this formula the operator T is the large gauge transformation 
operator which commutes with the Hamiltonian [T, H] = 0. The 
physical vacuum state is unique and constructed as a superposi- 
tion of \n) states. In QFT approach the presence of n different 
sectors in the system is reflected by summation over n £ Z in 
definition of the path integral in QCD. It should not be con- 
fused with conventional term "degeneracy" when two or more 
physically distinct states are present in the system. 



dispersive contribution, can be restored through the ab- 
sorptive part using conventional dispersion relations, and 
has a negative sign in accordance with general principles. 
This conventional physical contribution is saturated in 
this model by the lightest 77' field. It enters xqcd pre- 
cisely in such a way that the Ward Identity (WI) ex- 
pressed as XQcni'nT'q = 0) = is automatically satisfied 
as a result of cancellation between the two terms. If the 
contact non-dispersive term with "wrong sign" was not 
present in the system, the WI could not be satisfied as 
physical states always contribute with negative sign in 
eq. (6). 

One should note that the number of tunnelling events 
per unit time per unit volume (3) in pure gauge theory 
in this model (with no quarks) precisely concurs with 
the absolute value of the energy density of the system. 
Furthermore, the topological susceptibility in pure gauge 
theory calculated as the second derivative of Eym (^) with 
respect to precisely coincides with non-dispersive con- 
tact term with "wrong sign" explicitly and directly com- 
puted in (6). Indeed, 



-E'ym(^) 



Xym{0 = 0) = 



cos , 

L \Nc 



(7) 



d^EYu{0) 



96*2 



9=0 



NrL' 



where we keep only the lowest branch I — in expression 
for cos ( ^"^^' 1 to simplify formula (7), see detail discus- 
sions with complete set of references on this matter in 
[21]. In other words, the contact term in pure gauge the- 
ory xym = jYL '"^'^ ^^ interpreted in terms of number 
tunnelling events between different topological sectors in 
the system. Therefore, there is no surprise that it has 
the "wrong sign" as the relevant physics can not be de- 
scribed in terms of propagating physical degrees of free- 
dom, but rather, is described in terms of the tunnelling 
events between different (but physically equivalent) topo- 
logical sectors in the system. 



III. CASIMIR-TYPE BEHAVIOUR IN 
DEFORMED QCD 

Up to this point the theory was formulated on R"^ x S^ 
with small compactification size L for compact time co- 
ordinate S^ and infinitely large space K^ describing three 
other dimensions. As explained in section I, we are actu- 
ally interested in behaviour of the system when a space 
with large dimensions M'^ receives some small modifica- 
tions, for example the theory is defined in a ball R^^ -^ B"^ 
with L being a very large size of the compact dimension 
of the sphere S^ which is a boundary of the ball B'^ . Such 
a modification can be thought as a simplest way to model 
and test the sensitivity of our theory to arbitrary large 
distances such as size of our visible universe determined 
by the Hubble constant -ff/AqcD ^ lO"**^. We want 
to know how the topological susceptibility of the system 



which describes the 9 dependent portion of the vacuum 
energy i?vac(^ = 0) changes with slight variation of size of 
the system. We assume that L^i/^^^10 Gpc is much 
larger than any other scales of the problem. Essentially 
we want to see whether our deformed QCD model with a 
mass gap nia predicts an exponential scaling typical for 
a free massive particle 

A£:(L) = [E(M^) - E{R^)] - exp(-m^L) (8) 

or, it demonstrates a Casimir type behaviour 

AE{h)^[E{M^)-E{R^)]^^ + o(^^^ . (9) 

If we did not have a non-dispersive contribution in our 
system, we would immediately predict the behaviour (8) 
as the only available option for a gapped theory in close 
analogy with conventional Casimir computations for a 
massive particle AE{L,) ~ cxp(— mL), see e.g. review 
paper [17]. However, our system is much richer, more 
complicated, and more interesting, as it exhibits a non- 
dispersive term resulting from degeneracy of topological 
sectors in gauge theory as discussed in the text. This 
contact term, being unrelated to any physical degrees of 
freedom, may provide different scaling properties since 
conventional dispersion relations do not dictate its be- 
haviour at very large distances. As we shall argue be- 
low, the deformed QCD indeed exhibits the Casimir type 
behaviour (9) in drastic contrast with the conventional 
viewpoint represented by eq.(8). As we reviewed in sec- 
tion I we interpret a tiny deviation of the 9- dependent 
vacuum energy i?vac in expanding universe (in compari- 
son with Minkowski space-time) as a main source of the 
observed dark energy. The Casimir type behaviour (9) 
plays a key role in possibility of such an identification. 

We start our discussions in section HI A with con- 
ventional 4d instanton computations [22] in which in- 
frared regularization for some gauge modes is required 
and achieved by putting the system into a sphere with 
finite radius L. It allows us to compute power like cor- 
rections to the standard instanton density [22]. However, 
the corresponding corrections being computed for a fixed 
instanton size p can not be interpreted as a physically 
observable quantity because the integral / dp over large 
size instantons diverges for this system when semiclassi- 
cal approximation for large p breaks down. Nevertheless, 
this example explicitly shows when and why a Casimir 
type correction (to conventional formula computed in in- 
finite R^ space) emerges. 

Next, we compute a similar correction for the "de- 
formed QCD" model in section IIIB wherein a Casimir 
type correction also appears, resulting from the same 
physics related to topological sectors of the theory. In 
contrast with the previous case, the correction computed 
in this system is physically "observable" quantity as it 
represents the vacuum energy of the system. Indeed, 
the tunnelling transitions in this case are described by 
weakly coupled monopolcs, such that semiclassical com- 
putations of the vacuum energy (3), (5) expressed in terms 



of the density C of pseudo-particles are fully justified. The 
size of pseudo-particles (fractionally charged monopoles) 
which describe the tunnelling events in this model is fixed 
by construction see [20] [21] for the details. 



A. Casimir-type corrections for 4d instantons 

Our goal here is to study a power like correction to the 
instanton density described in the classic paper [22]. As 
such, we adopt 't Hooft's notation, and in particular, use 
the same background-dependent gauge C4 ~ Vfj^A^^^, 
which drastically simplifies all computations. Essentially, 
the problem is reduced to analysis of the normalization 
factors for finite number of zero modes (8 for SU{2) gauge 
group) in this gauge wherein the system is defined in a 
sphere with large but finite radius radius L. Essentially 
we follow the construction described in section XI of ref. 
[22]. The corresponding normalization factor explicitly 
enters the expression for the instanton density as it ac- 
companies the integration over collective variables. The 
contribution from non-zero modes does not exhibit such 
corrections, see the few comments on this issue at the end 
of this section. We now concentrate on the zero modes 
and power like corrections which accompany the normal- 
ization factors if the system is defined on a large but finite 
space Bl (four dimensional interior of a ball of radius L) 
rather than an infinite space M.'^. 

We start with four translational zero modes which have 
the form 

^r"M~'?'^M-(i+^')"'' ^ = i,-,4 (10) 

where we literally use 't Hooft's notations for rya^^ sym- 
bols and dimensionless coordinate r 



x^, measured in 



units of p = 1. Computing the corresponding correction 
factor due to the translation zero modes Ktr., we have 



Ktr 



/„'-d^x[A-q"(^)]^ 
/o°°d4^[Ar"(^)]^ 



3 

17 



^(17) 



.(11) 



The corresponding correction factor to the instanton den- 
sity has power like correction as anticipated. As a result 
of additional rotational symmetry one should expect, in 
general, L~^ corrections, while translation zero modes 
lead to a much smaller correction ~ L~"' as eq. (11) 
shows. It will be neglected in what follows. Dilaton and 
global gauge rotations lead to ^ L~^ as we discuss below. 
For the dilaton zero mode 



/la qu 



Va^.^x•' {1 + r"^)- 



a similar formula reads 



Kdii. 






-i 



^(17^ 



(12) 



,(13) 



such that the correction to the instanton density is pro- 



Computing the corresponding contribution due to 
three zero modes related to global gauge rotations re- 
quires much more refined analysis as explained in [22]. 
This is due to the specific features of the background de- 
pendent gauge C4 = P^Afj'^" when the corresponding 
three modes are pure gauge artifact. As shown in [22] 
the corresponding contribution is finite, but very sensi- 
tive to the infrared regularization determined by the size 
R of large sphere. The corresponding contribution to 
the instanton density is ^ (X^V)^^'^ where V is the four 
volume, while A4 ~ y^^ is defined as follows 



A4 = 



jyd^x[r{bw 



& = 1,2,3, 



(14) 



i'^'ib) = 'naf_cur]b^LX- 



(1+^2)' 

The corresponding power like corrections can be com- 
puted in a similar manner to the other zero modes, ex- 
cept that we must retain the regularization since the de- 
nominator above diverges as ^ V. So we have the two 
correction factors 



and 



l^de 



£d^x[r,{b)? 
j;^d^x[r,m 



^ viR)j^d''x[r{b)? 

V(l^) J^ d^x[-iP''{b)]^ 



1-^+^(17) 



4 1 

1-7^ + ^(77; 



L2 



X4 



The fraction, V{R)/V{h), is the correction to V in the 
instanton density factor, and is included here so that we 
can take the regularization i? — > cx). The combined gauge 
rotation correction factor is then 



'^rot. 



I^de 



^ + 17+^(17) 



(15) 



such that the correction to the instanton density is pro- 
portional to (Krot.)'^^^ ~ (1 + 2p')- Accidentally, for 
SU{2) gauge group the leading L"^ correction from the 
dilation (13) and global gauge rotations (15) exactly can- 
cel each other. This accidental cancellation does not hold 
for general SU{N) gauge group when power of K^ot. en- 
ters the instanton density with a different power. 

We remark here that the technique used in [22] is essen- 
tially a variational approach wherein the boundary con- 
ditions are implemented implicitly rather than explicitly. 
It allows us to use all the zero modes (10), (12), (14) as 
well as standard classical instanton solution in the origi- 
nal form defined on R* in which the conformal invariance 
is a symmetry of the system. So in this approach, neither 
the instanton itself, nor its zero modes (10), (12), (14) are 
solutions of the equation of motions which vanish at the 
boundary. This approach has been tested in many fol- 
low up papers, and we adopt it in the present work using 



the same technique in the next section. We also point 
out that the conformal invariance is exphcitly broken in 
the one instanton sector by the size of the instanton p, 

2 

such that corrections take the form (jj-)"- It is restored 
by the integration J dp. However, in this paper we are 
interested in by the computation in one instanton sector 
only when dimensional parameter p is explicitly present 
in the system, and it is small and fixed. 

The important message here is that such kind of power 
correction do appear in general. The source of these cor- 
rections is a long range tail of zero modes. We can not 
derive a definite conclusion from these computations be- 
cause the integral over large size instantons J dp diverges 
and the semiclassical approximation breaks down. How- 
ever, the same problem studied in the deformed QCD 
model considered in Section IHB does not suffer from 
such deficiencies as semiclassical computations are under 
complete theoretical control. Thus, a Casimir like correc- 
tion to the monopole fugacity ( in this model is explicitly 
translated to the correction to the vacuum energy density 
and topological susceptibility (7), supporting (9) and in 
huge contrast with naive expectation (8). It is impor- 
tant that the source of the corrections in the deformed 
QCD model is the same as in undeformed QCD consid- 
ered here, and that source is the long range tails of the 
zero modes, which lead to large distance sensitivity. The 
only difference is that the role of the instanton size p in 
computations above in the one instanton sector is played 
by the inverse monopole's mass m^ in next section IHB. 
Because it is a true scale of the problem however, tti^ is 
not integrated over as p is. 



B. Casimir- type corrections for 3d monopoles 

We now turn to the deformed gauge theory described 
in [20, 21] wherein the low-energy behaviour is given by 
a U{1)^ Coulomb gas of monopoles in Euclidean R^. 
Basically, we want to understand the dependence of the 
monopole fugacity, C, which comes out of the measure 
transformation to collective coordinates, on the size of 
the system, L. In this case, as in the previous section, 
we consider the interior of a sphere of large but finite ra- 
dius L. There are four zero-modes present in this system: 
three translations since the monopoles are in M.^, no di- 
lations since the monopole size is fixed by the symmetry- 
breaking scale in this model mw, and one gauge rotation 
since the gauge group for a given monopole is U{1). As 
in [22] , we work in a regular gauge to remain sensitive to 
the large distance physics. The monopole solution in the 
"hedgehog" regular gauge is given by 



''-'u^'''-' ~ ^fJ.i^a 2 



x\ 



1- 



mw \x\ 



sinh (mw \x\ 



(16) 



b^ix) = — 2 ["^w 1^1 coth {mw \x\) - 1] , (17) 



where we adapted notations from [23, 24] treating the 
monopole measure in supersymmetric Yang -Mills theory. 
In formula (16) w° denotes the three spacial gauge fields 
for the classical solution, and 0° the gauge field in the 
compact time direction (the "Higgs" field in this model) 
when all fields can be combined in a single 4d field «,„. 

We then want to compute the correction factors for 
the collective coordinate measure coming from these four 
zero modes when the system is defined in a large but finite 
sphere. We closely follow the 't Hooft's treatment [22] 
presented in previous section HI A. We start by consider- 
ing the translation modes defined by the spacial deriva- 
tive of the classical monopole solution (16) 



Z^^{l^) = -dXnix -Z)+ V„,vl 



(18) 



where second term on the right hand side of eq. (18) is 
necessary to keep Z'^^{v) in the background gauge, see 
[23, 24] for the details. This leads us to the following 
expression for correction factor due to the translation 
zero modes 



Ktr. 



/>^x[Z°,(^)]^ 
j^ d^x[Z'^^M)] 



1 



' +0(A) 



mw^ 



^L2 



(19) 



Next we consider the gauge rotation zero-mode. As 
in the previous section, the contribution to the collective 
coordinate measure, and so the monopole fugacity, is ~ 
(Ay) 2 where V is the three- volume and A is given by 



A 



!vd'x[Bl\' 
/^d3a;[0'^]2 



(20) 



Bl = -e^.pa.< 



V„ 



Again, the denominator diverges as ~ y and we look at 
the two correction factors 



and 



'*do 



S^d^^[B^,Y 



V{R)j,d^xm 



1 



mv^L 



-^(17) 



niw^ 



1 



o{,,. 



The total correction factor for the gauge rotation mode 
is then 



'>^rot. 



^dc 



1 



mw^ 



-^^ 



(21) 



and therefore the total correction to the monopole fugac- 
ity from the (20) is ^/K^. ~ (1 + |^). Assembling the 
total correction to the fugacity. 



3/2 1/2 



1 



1 



2mw^ 



^(17^ 



(22) 



Thus, the deformed QCD, when put on a manifold with 
a boundary, receives some corrections to the monopole fu- 
gacity compared to Minkowski space that are power-like 



in the manifold size. The correction (22) to the monopole 
fugacity leads immediately to the same correction to the 
topological susceptibility and so the background energy 
density through the relation (7). To be more precise, 



C(IL) = C ■ 



^ 2mwL^^^L2^ 



(23) 



where ^ is the monopole fugacity which enters the rela- 
tion (7) computed in infinite Minkowski space. We em- 
phasize that the energy density changes in the bulk of 
space-time, not only in the vicinity of the boundaries, 
similar to the Casimir effect when the bulk energy den- 
sity changes as a result of merely presence of the bound- 
ary. To reiterate, the deformed QCD, despite the pres- 
ence of a mass gap, displays a suprising Casimir-like sen- 
sitivity to large distance boundaries, such that the en- 
ergy density differs from the Minkowski space value by 
/S.E ^ — ^r. Again, this is in contrast to the naive ex- 
pectation based on analysing the physical degrees of free- 
dom, /S.E ^ g-^™L with m ~ rria being the lowest mass 
scale of the problem (8). 



C. A few general comments 

Computations of the Casimir corrections presented 
above were based on an analysis of the zero modes when 
corresponding normalization factor explicitly enters the 
instanton/monopole density. Now, we want to present 
some arguments suggesting that a corrections due to 
the non-zero mode contributions can be neglected, and, 
therefore, it cannot cancel the zero modes contribution. 
Indeed, the computation of non-zero mode contribution 
is reduced to analysis of the phase shifts in the scatter- 
ing matrix which can not change the normalization of 
the wave function itself as the only changes occur are 
the phase shifts. An absolute normalization is dropped 
from the final formula for the instanton/monopole den- 
sity when the ratio of the eigenvalues is considered. This 
argument is consistent with observation that non-zero 
mode contribution depends on matter context of the the- 
ory as it varies when massive scalar of spinor fields in dif- 
ferent representations are part of the consideration. At 
the same time, the Casimir type corrections computed 
above are exclusively due to the gauge portion of the 
theory, not its matter context. Indeed, these Casimir 
corrections were derived in pure gluodynamics. So, it is 
difficult to imagine how a Casimir correction to non-zero 
mode contribution (if it is nonzero) may cancel a Casimir 
type correction originated from analysis of gauge zero 
modes. 

We also comment that the correction L~^ occurs as a 
manifestation of a slow power like decay of the zero modes 
in the background of a topologically nontrivial gauge con- 
figuration. It should be contrasted with conventional be- 
haviour of zero modes with a mass gap present in the 
system from the very beginning (for example, the well 
studied problem of a double well potential). In former 



case, the zero modes decay according to the power law 
and leads to the Casimir type correction, while in the 
later case, the zero modes are well localized configura- 
tions which decay exponentially fast at large distances 
and can not be sensitive to large distance physics. The 
mass gap is present for all physical degrees of freedom in 
both models. However, in the former case the mass gap 
emerges as a result of the same instanton/monopole dy- 
namics, while in the later a mass gap was present in the 
system from the very beginning and it was not associated 
with any instanton/monopole dynamics. QCD obviously 
belongs to the former case, and we therefore expect this 
effect will persist in real strongly coupled QCD. 

Next, our computations of the Casimir correction to 
the instanton/monopole density are based on assumption 
of the dilute gas approximation. This is enforced in sec- 
tion III A by a finite instanton size p which is kept fixed 
and small. On other hand, the semiclassical approxima- 
tion in section III B is automatically justified due to para- 
mctrically small fugacity C, and total neutrality in this 
system is automatically achieved as long as the size of the 
system L is much larger than the Debye screening length 
m^^, see (4). In other words, we assume L ^ m~^ such 
that neutrality of the system is automatically satisfied 
with exponential accuracy. The finite size of the manifold 
does not spoil this neutrality if condition L ^ rn"^ is sat- 
isfied. Furthermore, the computation of the monopole's 
fugacity Q and corresponding corrections (23) can be per- 
formed without taking into account of the interaction of 
a monopole with other particles from the system as it 
would correspond to higher order corrections in density 
expansion ~ ^2. This is precisely the procedure which 
was followed in the original computations by Polyakov in 
[25] and in the deformed QCD model in ref.[20] at weak 
coupling. 

Also, we emphasize that in variational approach devel- 
oped in [22] neither the classical solution nor the corre- 
sponding zero modes vanish at the boundary of a finite 
size manifold. The constraints related to the finite size L 
of the manifolds are accounted for implicitly rather than 
explicitly in this approach. In particular, one should not 
explicitly cut off the classical action of the configuration 
as a result of finite size L where instanton/monopole is 
defined as this contribution is implicitly taken into ac- 
count by variational approach. However, even if we use 
an explicit cutoff for classical solution (as some people 
suggest) it still cannot cancel the zero mode corrections 
as these terms have different behaviour in N . Indeed, 
the correction to the classical solution would be one and 
the same for any N ^ while corrections due to zero modes 
depend on TV as correction (21) counts number of gauge 
rotations for SU{N) gauge theory. 

Finally, it is quite possible that we overlooked some 
other possible corrections (for example, some corrections 
due to the boundaries which may occur in close vicinity 
of these boundaries). We emphasize that our main re- 
sult is not a computation of a specific coefficient in front 
of the correction to fugacity in equation (23). Rather, 



our main point is that these type of corrections do oc- 
cur in a system with a gap, and it is very difBcuh to 
imagine that some boundary corrections may mysteri- 
ously cancel these computed bulk corrections (as some 
people suggest). Therefore, we present below some ar- 
guments and examples suggesting that Casimir type be- 
haviour in gauge theories is in fact quite generic, rather 
than a peculiar feature of the system. 



IV. TOPOLOGICAL SECTORS AND THE 
CASIMIR CORRECTION IN QCD 

In this section we want to present few generic argu- 
ments suggesting that the emergence of a Casimir-like 
behaviour is not an accident, and not a computational 
error. Rather, the effect has a deep theoretical roots as 
argued in [26] . We review these arguments starting with 
analogy with the well known Aharonov-Casher effect as 
formulated in [27]. The relevant part of that work can be 
stated as follows. If one inserts an external charge into 
superconductor when the electric field is exponentially 
suppressed ^ exp(— r/A) with A being the penetration 
depth, a neutral magnetic fluxon will be still sensitive to 
an inserted external charge at arbitrary large distance. 
The effect is purely topological and non-local in nature. 
The crucial point is that this phenomenon occurs, in spite 
of the fact that the system is gapped, due to the presence 
of different topological states in the system. We do not 
have a luxury of solving a similar problem in strongly 
coupled four dimensional QCD analytically. However, 
one can argue that the role of the "modular operator" 
of [27] (which is the key element in the demonstration of 
long range order) is played by large gauge transforma- 
tion operator T in QCD which also commutes with the 
Hamiltonian [T, H] = 0, such that our system must be 
transparent to topologically nontrivial pure gauge config- 
urations, similar to transparency of the superconductor 
to the "modular electric field" , see [26] for the details. 

We interpret the computational results in a number of 
systems where Casimir like corrections have been estab- 
lished as a manifestation of the same physics which can be 
described in terms of the operator T. We should mention 
that there are a few other systems, such as topological 
insulators, where a topological long range order emerges 
in spite of the presence of a gap in the system. We refer 
to ref [26] for relevant references and details. 

There arc a number of simple systems in which the 
Casimir type behaviour AE ^ L^^ + 0{h)~'^ has been 
explicitly computed. In all known cases this behaviour 
emerges from non-dispersive contributions when the dis- 
persion relations do not dictate the scaling properties of 
this term. 

The first example is an explicit computation [28] in 
exactly solvable two-dimensional QED defined in a box 
size L. The model has all elements crucial for present 
work: non-dispersive contact term which emerges due 
to the topological sectors of the theory. This model 



is known to be a theory of a single physical massive 



field. Still, one can explicitly compute AE 



L- 



contrast with naively expected exponential suppression, 
AE ~ e~^ [28] . Another piece of support for a power like 
behaviour is an explicit computation in a simple case of 
Rindler space-time in four dimensional QCD [3, 18, 29] 
where Casimir like correction have been computed us- 
ing the unphysical Veneziano ghost which effectively de- 
scribes the dynamics of the topological sectors and the 
contact term when the background is slightly modified. 
Thus, power-like behaviour is not a specific feature of 
two dimensional physics as some people (incorrectly) in- 
terpret the results of [28]. 

Our next example is 2d CP^~^ model formulated on 
finite interval with size L [30]. In this case one can ex- 
plicitly see emergence of AE ^ L^^ in large N limit 
in close analogy to our case (23) where a theory has a 
gap, but nevertheless, exhibits the power like corrections. 
The correction computed in [30] also comes from a non- 
dispersive contribution which can not be associated with 
any physical propagating degrees of freedom, similar to 
our case (23). 

Power like behaviour AE ^ L^^ is also supported by 
recent lattice results [31]. The approach advocated in 
ref. [31] is based on physical Coulomb gauge, in which 
nontrivial topological structure of the gauge fields is 
represented by the so-called Gribov copies leading to a 
strong infrared singularity. Thus, the same Casimir-like 
scaling emerges in a different framework where the un- 
physical Veneziano ghost (used in refs. [3, 18, 29]) is not 
even mentioned. 

The very same conclusion also follows from the holo- 
graphic description of the contact term presented in [26] . 
The key element for this conclusion follows from the fact 
that the contact term in holographic description is de- 
termined by massless Ramond-Ramond (RR) gauge field 
defined in the bulk of 5-dimensional space. Therefore, it 
is quite natural to expect that massless R-R field in holo- 
graphic description leads to power like corrections when 
the background is slightly modified. 

To avoid any confusion with terminology we follow [26] 
and call this effect as "Topological Casimir Effect" where 
no massless degrees of freedom are present in the sys- 
tem, but nevertheless, the system itself is sensitive to 
arbitrary large distances. It is very different from con- 
ventional Casimir effect where physical massless physical 
photons arc responsible for power like behaviour. From 
the holographic viewpoint discussed in [26] the "Topo- 
logical Casimir Effect" in our physical space-time can 
be thought as conventional Casimir effect in multidimen- 
sional space when massless propagating R-R field in the 
bulk is responsible for this type of behaviour, although 
this field is not a physical asymptotic state in our four 
dimensional world. 



V. CONCLUSION AND FUTURE DIRECTIONS 

We tested a sensitivity of a deformed QCD model 
with non-trivial topological features to arbitrary large 
distances. A naive expectation based on dispersion re- 
lations dictates that a sensitivity to very large distances 
must be exponentially suppressed (8) when the mass gap 
is present in the system. However, we argued that along 
with conventional dispersive contribution there exists a 
non-dispersive contribution, not related to any physi- 
cal propagating degrees of freedom. This non-dispersive 
(contact) term with the "wrong sign" emerges as a result 
of topologically nontrivial sectors, and can be explicitly 
computed in our model. The variation of this contact 
term with variation of the background leads to a power 
like "Topological Casimir Effect" (9) in accordance with 
other arguments presented in section IV and in contrast 
with the naively expected exponential suppression (8). 

The "Topological Casimir Effect" in QCD, if confirmed 
by future analytical and numerical studies, may have pro- 
found consequences for understanding of the expanding 
FLRW universe we live in. We already mentioned in sec- 
tion I that the observed DE (1) may is fact be just a 
manifestation of this "Topological Casimir Effect" with- 
out adjusting any parameters. In the adiabatic approxi- 
mation the universe expansion can be modeled as a slow 
process in which the size of the system adiabatically de- 
pends on time L(i) which leads to extra energy as equa- 
tions (9, 23) suggest. Such a model is obviously con- 
sistent with observations if L(t) is sufficiently large [32]. 
We do not insist that this is the model of our universe. 
Rather, we claim that if the effect persists in strongly 
coupled QCD, the energy density which can not be iden- 
tified with any physical propagating degrees of freedom, 
is sensitive to arbitrary large distances as a result of non- 



trivial topological features of QCD. Different geometries 
(such as FLRW universe) obviously would lead to dif- 
ferent coefficients. Nonetheless, the important message 
from these computations in our simplified model is that 
the energy density in the bulk is sensitive to arbitrary 
large distances comparable with the visible size of the 
universe, and that this sensitivity comes not from any 
new physics but simply from the proper treatment of the 
topological structure of QCD. 

We add that a comprehensive phenomenological anal- 
ysis based on this idea has been recently performed in 
[33] where comparison with current observational data 
including Snia, BAO, CMB, BEN has been presented, see 
also [29, 34-38] with related discussions. The conclusion 
was that the model (1) is consistent with all presently 
available data, and we refer the reader to these papers 
on analysis of the observational data. 

Finally, what is perhaps more remarkable is the fact 
that the "Topological Casimir Effect" which is the sub- 
ject of this work can be, in principle, experimentally 
tested in heavy ion collisions, where a similar environ- 
ment can be achieved, see [18, 19] for the details. In 
particular, the V odd correlations observed at RHIC and 
LHC have been interpreted in [18, 19] as a result of long 
range order represented by the "Topological Casimir Ef- 
fect". 
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